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&*• Abstract: 



It is proven that the energy of a quantum mechanical harmonic oscillator with 



a generically time-dependent but cyclic frequency, u>o{to) = ^o(O), cannot decrease 
on the average if the system is originally in a stationary state, after the system goes 



through a full cycle. The energy exchange always takes place in the direction from 



the macroscopic system (environment) to the quantum microscopic system. 
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1. A Theorem 

Many physical systems reduce effectively to that of a harmonic oscillator with a 
time-dependent frequency (we set m — 1), 

H(t) = £ + ±u>(t)W, (1-1) 

which is thus of certain interest [1]. In particular, the case of a periodic (cyclic) 
variation of uo such that 

u(t )=u(p) (1.2) 

is of great interest. We shall not specify the time variation oiu(t) otherwise. Suppose 
furthermore that the system is initially in one of the stationary states, 1^), with energy 
E in . We prove below, independently of how the system goes through the cycle, that 

E fin = E(t ) = (iJ>(t )\H\i!>(t )) > E in , (1.3) 

i.e. on the average the microscopic system can never give excess energy to the external 
environment. Such a result might appear surprising at first sight, since during a 
generic time variation of external parameters, the system can either give the energy 
away to or absorb it from the environment. Also, such a theorem certainly does not 
hold in general in a system with a finite number of independent states, such as a spin 
| system in a varying magnetic field. 

The proof is easiest in the Heisenberg picture. Heisenberg equations of motion 
tell us that q(t), p(t) are linear combinations of q, p; they furthermore preserve the 
commutation relation (we set h — 1) 

[q(t),p(t)]=i. (1.4) 

In other words the time evolution is described by an Sp(2) transformation 

y(t) a = S a p(t) y(0)/3, S aj3 S l5 £p 5 ^e ai SeS T = e, (1.5) 

where we wrote y a = (q,p). Or writing 



S — J r ) , a,b,c,d real 
the condition (1.5) reduces to 




ad-bc = det(S) = l. (1.6) 



The energy expectation value at time t is given by 

E(t) = (mm) = l -^wf+u{tf q {tf\i,). (i.7) 

At the end of a cycle (we set u;(0) = 1), 

Efin = l^\PH(t ) 2 + q H (to) 2 \ip)- (1.8) 



By using Eq. (1.5) and by defining 



1 
D = ^(iP + PQ) 



the quadratic form becomes 



(S u q + S 12 p) 2 + (S 2 iq + S 22 p) 2 = 

= (Su + S 2 21 )q 2 + (5 2 2 + S 2 22 )p 2 + 2(5-11512 + 5 21 5 22 ) J D. 

By using the Virial theorem 

\(i>\p 2 \i>) = \m 2 H) = \e iu (1.9) 

and the fact (valid for real wave functions) that 

ty\D\f)=0 (1.10) 

the final energy is given by 

Efin = -E in [Su + S 21 + 5^2 + S 22 ] . (1-11) 

The problem is then to find the minimum of a quadratic form Q = a 2 + b 2 + c 2 + d 2 
under a constraint ad — be = 1. Upon introduction of a Lagrange multiplier, the 
extremum of a 2 + b 2 + c 2 + d 2 + 2\(ad — be — 1) is given by 1 

A = — 1, a = d, b = — c, 5=1 ) , a 2 + 6 2 = 1, 

that is, when 5 is orthogonal. In that case the quadratic form takes the value 2 and 
therefore in general 

Efin = - E in Q > E in . (1.12) 



-'^As the set of the evolution matrices is unbounded, this extremum can only be a minimum. In 
an alternative proof given in (C.2) this fact is obvious. 



Remarks 

(i) In the adiabatic limit, the system "follows" the variation of the spectrum while 
staying in the "initial" eigenstate, and comes back to the original state, so we 
expect Ef —* E in . In the sudden limit, the wave function does not make it to 
change as the external parameter goes through a (too) quick cyclic variation, 
so that Ef — ► E in again. 

(ii) The above result is valid also in the case of a forced oscillator. Consider 

H=^(p 2 + uj(t) 2 q 2 )- K (t)q, (1.13) 

where n(t) is an arbitrary function with «(0) = k(T) = 0. The Heisenberg 
equations are 

p = i[H,p] = —u 2 q + k 

q = t[H,q}=p (1.14) 

so 

q = —uj 2 q + k. (1.15) 

Let us now consider a solution Q c (t) of (1.15) with the boundary condition 

Qe(0) = 4(0) = 0; 

it follows immediately that q is a sum of the homogeneous solution plus the 
particular one Q c 

q H (t) = Q c (t) + q ( r°\t), 

p H (t) = Qc(t)+p ( H = °\t). (1.16) 

One has then 

E f = 4 K=0) + \ [Q 2 + Q 2 ) + (Qc(q { r°\t)) + Qc{p { H = °\t))). 

As the expectation values of qu and pn vanish in the initial stationary state it 
follows that 

Ef = Ef= Q) + l -(Q 2 + Q 2 )>E l , (1.17) 

where the result of the preceding paragraphs has been used. 



(iii) In perturbation theory the theorem can be easily seen to be valid. By writ- 
ing u(t) 2 x 2 = u 2 x 2 + 5u(t)x 2 , 5u(t ) = 5u(0) = 0, the first-order transition 
probability 



p "-¥ 



z 



(1.18) 

is larger for the process n — »■ n + 2 than for n — > n — 2, hence (H) > E in . 



(iv) In perturbation theory, one can actually verify the theorem for a more gen- 
eral class of perturbing potentials. Indeed, it can be easily seen that for any 
perturbation of the form, 

AV(t, x) = 8u(t) x N , 8u(t ) = 5lo(0) = 0, (1.19) 

the theorem is valid, as \(x N ) n+mjn \ > \(x N ) n - m ,n\ (m > 0). 

This would suggest that actually the theorem holds for a wider class of periodic 
potentials, V(t ,x) = V(0,x). 

(iv) The theorem obviously applies for a system of iV independent oscillators 



N 2 1 
H=Y\ P -k + ^{t) 2 q 2 l (1.20) 



^ L 2 2 

i=l 

with arbitrary, periodic variations of u, Wi(t ) = c<->j(0), if the initial state is in 
a stationary state |ni, ri2, . . . , n^). 

(v) The theorem cannot hold for a generic initial pure state 2 . Since in some cases 
the final energy expectation value (i.e., in the pure state U(to) |^(0))) is strictly 
higher than the initial value (Sec. 2.), it suffices to consider the time- reversed 
process of such an evolution, to find a counter example 3 . 

On the other hand, for an initial state which is not an eigenstate of energy, the 
operational meaning of the theorem itself would become somewhat unclear. 

(vi) There is an important case in which the theorem applies for a initial mixed 
state. Consider the N oscillator system of (1.20) and suppose that the system 
is originally in the thermal equilibrium with a heat bath at temperature T. The 
theorem then clearly applies in a statistical and quantum average sense. 



2 We have however found a sufficient condition for the theorem to be valid: for initial states of 
the form -0 = J^ n a n il> ni n = uq + Am, m € Z, with a fixed no, the theorem can be shown to hold by 
a slight generalization of the proof given here. 

3 We thank Tomas Tyc for a useful communication on this point. 



(vii) The energy gain factor, R = {if)\H\ip)/E in , is universal, in the sense that it does 
not depend on the particular initial stationary state chosen. 

(viii) Note that the (classical) canonical equations of motion and Heisenberg equa- 
tions have the same form, and the evolution matrix S are the same in both 
cases. As a result, an analogous theorem holds in classical mechanics, if one 
takes an average over random initial values (p, q) over a given (fixed-energy) 
trajectory. 

2. Example: Inverse Linear Variation of the Frequency 

Consider the oscillator Eq.(l.l) with frequency varying as 

w (*) = ^» \(t) = l + vt, (2.1) 

where uq adn v are constants. From Heisenberg equations one gets 

By setting q(t) = A(t) /3 , one gets a characteristic equation 

p (P - 1) + n 2 = o, 

with solutions 

fa = l±8, 5= \J\-V\ (2.2) 

so that the general solution of the Heisenberg equations reads 

q H = Cl \^ + c 2 \P\ p H = v (ciftA* + c 2/ 3 2 A A ) . (2.3) 

By imposing the initial condition one finds 

qn{t) = J^j(ftAA-ftA*) g +i(A*-Afr)p}, 

p H (t) = ^ U fr (A ft - A ft ) q + I (ft A' 32 - A A ft ) p} . (2.4) 

As a check, consider the adiabatic limit, v — > 0, Q — ► oo. One has /^ — iicJo/^j so 
that 

A^m _> (l + t ,t)^/^ e T !W » f , 
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and 

quit) -> - {e~ iW{)t + e iuJot ) + -^—(-e' iuJot + e iuJot ) = qcoscu t + — smu t 
2 2icu too 

which is the correct result. 

Writing Eqs.(2.4) in the form of Eq.(1.5) with S = S(uo, v, A) one gets, by inserting 
this in Eq.(1.7) and by using the Virial theorem, 



W> = \e* 



\ 2 Wll + ^12) + ^21 + ^22 
,25 1 \-28 




E \ 2d + \' 2d + 2(A5 2 - 1) 
" T 4A^ • [ b) 

At time t, A = 1 + v t, the energy mean value (2.5) can be larger or smaller than 
the original energy depending on the sign of the velocity v (hence whether the scale 
factor A is smaller or larger than unity). 

However, as we are most interested in cyclic variations of 00, let us consider the 
the evolution from the original frequency uj to a final frequency ^, and then back 
to ujq. The second half of the evolution is simply described by the transformation 
S = S(^f-, —v, -r) so that the total evolution is described by the Heisenberg evolution 

= ^ c I qH f^ ) , ^ c = S(^, -v, 1) • S(u , v, A) (2.6) 

V Ph(0) J A A 

where S is defined by Eq.(2.4). 

The results for more general frequency variations are given in Appendix A and in 
Appendix B. 

In the case of the linear variation (2.5), (2.6), we have analysed numerically the 
energy gain factor R(u ,v, A) = (ip\H\ip) / E — | Tr [S cyc (S cyc ) T ] for various values 
of uj , v and A, A being the scale factor (iv = ujo/X) at the moment of the maximal 
contraction (A < 1) or expansion (A > 1). We find R(ujq, v, A) > 1 always as expected, 
but find also that: 

(i) At fixed A, R(ujq,v, A) — ► 1 both in the adiabatic (v — ► 0) and impulse approxi- 
mation (v — >• 00) limits, as expected; 

(ii) There are various resonance effects at small v (see Fig. 1); in particular, as a 
function of the velocity v at a fixed A, R reaches a maximum of order of 0(1/ X) 
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(for A < 1) or 0(A) (for A > 1), and then rather smoothly approaches the 
impulse-approximation value 1 asymptotically; 

(iii) When the system goes through iV cycles, the maximal energy gain factor (which 
occurs at certain critical velocity) behaves as R oc {\/X) N or R oc A^, a huge 
factor if A is large; 

(iv) There are values of (v, A) at which R attains values either exactly equal to or 
very close to unity (Fig. 1). 

(v) As a function of A, the maximum of R grows indefinitely as A — ► or as A — ► oo. 





Figure 1: Energy gain factor as a function of v for fixed A(= 10), at small and large v. 



3. Planck Distribution inside an Oscillating Cavity 

As a possible physical application of our considerations, let us consider the black 
body spectrum confined in a perfectly reflecting three-dimensional box of linear di- 
mension L. The electromagnetic fields are described by the Lagrangian 



where k n = irn/L. A redefinition of Qi n and a Legendre transformation leads to the 
Hamiltonian per unit volume 



ff = E(TP» + k2 *) + ^(T p ? 2 » + k2q: 



At temperature T the energy distribution is 



87r hv 2 

Suppose that at certain moment the box starts to contract or expand with a 
constant linear velocity. How does the energy distribution change with the linear size 
of the box? 

At time t — various modes are distributed according to the Planck distribution, 
(3.1). Each mode simply transforms as in Eq.(2.4): (\S\ ^$> 1, S — i\S\), 

E . H) _E n X a + \-* + 2(4#-l) ^E n 



Since 



u = Zjp£ ? L(t) = Lo(l + vt), \V\ = L \v\ <c, 



„ ijj n n c . /I ^ n 7r n c 

fi = - = ^»^ 5 = Vi- fi " 2 m (3 - 3) 

the process is adiabatic for all modes. The distribution remains Planckian; each mode 
is either red-shifted (in an expanding box) or shifted towards ultraviolet (a contracting 
box), simply by the scale factor A. 

Actually, the system is not so interesting as an application of the general theorem 
mentioned in Section 1, as it is entirely adiabatic; the only reason we discuss it here 
is its possible relevance to the so-called sonoluminescence phenomenon. 

In the single-bubble sonoluminescence [2], a gas bubble trapped in a liquid is 
made oscillating radially by acoustic waves, and at the moments its radius attains 
the minimum the bubble emits pulse of light in the visible to ultraviolet wavelength 
range. In a typical experiment, the linear dimension of the bubble contracts up to 
a factor of A = 10~ 4 : this is compatible with the observed increase of an effective 
temperature of a factor 10 4 or more. The energy density increases according to the 
Stefan's law, 

U = aT 4 ; a = 7.64 • l(r 15 ergcnr 3 K- 4 (3.4) 

but since the volume itself decreases by a factor A 3 the total energy of blackbody 
radiation increases only by a factor j-. It is possible that this excess energy is released 
in the form of visible light pulse at each cycle 4 . 

4 In an experiment at Laurence Livermore Ntional Laboratory [2] 10 5 — 10 6 photons are observed 



4. Conclusion 

A harmonic oscillator with a time-dependent frequency (or an ensemble of such 
oscillators) experiences a full cycle, with the system originally in a stationary state. 
After the cycle, the energy expectation value is predicted never to decrease, inde- 
pendently of the way the parameters vary with time during the cycle. Our result 
is somewhat reminiscent of the law of entropy increase, but concerns the evolution 
of quantum mechanical pure states, and no information loss is implied. A generic 
non-adiabatic "disturbance" by the external force always does work on the micro- 
scopic system on the average, increasing its energy. The energy flow is always in 
the direction from the macroscopic system to the microscopic one. In other words, 
the quantum mechanical harmonic oscillator cannot act as a perpetual machine, nor 
produce a net energy gain. A spaceship cannot continue her journey forever, getting 
its energy supply from the inexhaustible zero-point energy of the vacuum [3]. 
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Appendix A General Power Dependent Frequencies 



Consider a case of a generic power-behaved frequency, 

z=l + vt. 



V ; 2m 2 ° 



(A.l) 



For k = we recover the case discussed in the text. In the following, we shall set 
m — 1,Uq — 1. Note that v has the dimension of a frequency, so that to recover the 
dependence on u it suffices to replace it by v/ujq. The Heisenberg equation of motion 
gives 



Qh + z k 2 q H = 0. 



Multiplying this by 1/v 2 one gets 

dz 2 
Its general solution is of the form, 



— ?fl + ~^z k 2 q H = 0. 



2i k 2i k 

q H (t) = A^/zJ l/k (—z2) + By/zY 1/k (—z*] 
kv kv 



(A.2) 



(A.3) 



where J, Y are the Bessel functions of the first and second kind, respectively. Differ- 
entiation with respect to time yields 

Pn(t) = v J^<lH(t). 
The coefficients A and B are determined by imposing the initial conditions 

qu\ z=1 =q Ph\ z=1 =P- 



kv 
kv 



q-{Y 1+l {-) + {p-vq)-Y,{- 



q '^+\^ + ( p - V( ti' J ^— 



2 

k^kv 

2 
k y kv 



The expression for 

is quite complicated but can be analyzed numerically. We note that: 



E(t) = (0\±p H (t) 2 + l -z k - 2 q H (t) 2 \0). 



(A.4) 
(A.5) 
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1) The case k = —2 reduces to the known result; 



2) The behavior is qualitatively similar for all k, with a phase of monotonous 
increase following the initial, oscillating phase, as a function of v (see Fig. 2); 

3) At fixed final "physical scale", Zf in = A - *^, z k ,~ 2 = A -2 , the asymptotic 
behavior in v turns out to be always the same as in the case k = —2, i.e., 

E f ^^(l + y 2 ). (A.6) 

This assertion is based on a graphical evidence (see Fig. 2), for the moment 5 . 
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Figure 2: Energy as a function of V for a fixed A = 1/10, with k = —2, —3, —4. Expected asymptotic 
value is 1/4(1 + 100) = 25.25. 



Appendix B Exponential Dependence 



As a second example, let us consider the Hamiltonian, 

H =-p 2 + -e 2vt q 2 . 
T 2 H 



(b.i; 



The Heisenberg equation is 



q H + e 2vt q H = 0. 



5 For generic order, the functions J, Y have essential singularities at v — oo, which prevent us 
from analyzing either adiabatic or large v limit analytically. 
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Changing the variable to r = v t one has 

d I 



J.vt. 



which has 



q H (t) 
p H (t) 



-r^qu + -^e zvt q H = 
dr z v z 



AJ (-e vt )+BY (-e vt ) 

v v 

d ^L = -Ae^M-e 01 ) -Be vt Y 1 (-r 
dt y v ' y v 



vt\ 



(B.2) 



as general solutions. As before the coefficients are determined by the initial condition: 



2v 



pY (-) + qY 1 (-) 

v v 



B=^ 

2v 



p JoH + q Ji{-) 

V V 



The mean energy can be computed as in the preceding cases. 
E(t) 



vrV 



, - 2 j (-)y (-; 

107T Z V V 
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1 


Jo 


(-)Y (-) + J 1 (-)Y 1 

V V V 


V 



+Jl 



■v' 



<)+^: 



+ Ji> 



v 



Y 2 



1. 



V 



Yt 



+2Ji(-)^(-; 

v v 



M Vo(-) + M Vi(- 

V V V V 



+ 



+ 



J, 



°tf +J ^' 



Y 2 



,Z S 



Yl 



,z s 



(B.3) 



where z — e . It is also possible to get the asymptotic behavior in v at fixed z: 



E - = k< 1+ *-ib<-' " 



(, 2 - 1)= 



(B.4) 



which is compatible with (A. 6) asz = 1/A. 



Appendix C Creation and Annihilation Operators 



The whole problem can be analyzed by use of creation and annihilation operators. 
We introduce at each instance the variables qi(t), Pi(t) in which the frequency is 
diagonal, u>ij(t) = SijUi(t); then define di(t), a\(t) in the standard manner in terms of 
qi(t), Pi(t). The time evolution introduces a linear transformation among aj(£),aj(£), 
which has a general form, 



Aika k + Bikdk, 



a, 



A* k a k + B* k a\. 
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The coefficients must be such that the canonical commutation relations are preserved 
(in a matrix form): 

AA*-BB* = 1,AB T -BA T = 0. (C.l) 

For a single oscillator, the theorem of Section 1. can be immediately proven: 

E f = -+u{n\(a') ] a'\n) = -+uj{n\(\A\ 2 + \B\ 2 )(Ja\n) = 

= ^ +n (l + 2\B\ 2 )u>uj(n+^) = E m . (C.2) 

For iV oscillators which are originally in a stationary state (ignoring the zero point 
energy), 

Ei = (\P| y^^qto^) = y^^rij 

£, = (*| ^^a'taJl*) = (*| Y^WiiAikal + B ik a k )(A* eai + B* t a\)\^>) = 

i i 

= E ^ (A ik A* k n k + B ik B* k (n k + 1 ) ) . 

i . k 

For the diagonal matrices A e B (uncoupled oscillators), (C.2) trivially generalizes. 
For coupled oscillators, the generalization is not obvious. 

We content ourselves here with a somewhat weaker result: it states that the 
occupation number (number of phonons) does not diminish under cyclic variations of 
frequency. Indeed, 

N i = (*lE 2 0i l*} = J2 ni 

i i 

N f = (*| XyjoW = <*l £(4*4 + B lk a k )(A* e a e + 5>{)|*> = 

i i 

= J2( A ^ A >x + B ikB* k (n k + 1)) 

■i . k 

Eq.(C.l) implies 

E A ^ = E B * s * + 1 

i j 

therefore 

k i,k 
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